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Abstract 

The initial structure and axial decay of an array of 
streamwise vortices embedded in a turbulent pipe bound- 
ary layer is experimentally investigated. The vortices are 
shed in counter-rotating fashion from an array of equally- 
spaced symmetric airfoil vortex generators. Vortex struc- 
ture is quantified in terms of crossplane circulation and 
peak streamwise vorticity. Flow conditions are subsonic 
and incompressible. The focus of this study is on the ef- 
fect of the initial spacing between the parent vortex gen- 
erators. Arrays with vortex generators spaced at 15 and 
30 degrees apart are considered. When the spacing be- 
tween vortex generators is decreased the circulation and 
peak vorticity of the shed vortices increases. Analysis 
indicates this strengthening results from regions of fluid 
acceleration in the vicinity of the vortex generator array. 
Decreased spacing between the constituent vortices also 
produces increased rates of circulation and peak vorticity 
decay. 

Introduction 

The use of surface-mounted vortex generators is 
now standard practice in a variety of fluids engineering 
applications where a certain amount of flow conditioning 
is critical for the proper performance of flow components. 
Vortex generators are, at the most basic level, a protrusion 
into a fluid stream that sheds a trailing vortex (or vortices) 
into the boundary layer downstream. The key to the 
effectiveness of the vortex generator lies in the convective 
nature of this shed vortex. Two important mechanisms 
of “convectiveness” are identified: 

1. Mixing. The swirling motion of the fluid in the vor- 
tex mixes high energy fluid of the frees tr earn or core 
flow into the slower moving fluid of the boundary 
layer. Thus mixing originates with the normal (up 
and down) component of the crossflow induced by 
the vortex. In the downflow region of the vortex 
the boundary layer is thinned and local values of 
skin friction are correspondingly high. In the upflow 
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region, the boundary layer is thickened and local 
values of skin friction are correspondingly lowered. 
The overall or integrated effect of vortex mixing is to 
promote re-energization of the boundary layer fluid 
and boundary layer attachment. The mixing mecha- 
nism of vortex generator effectiveness is long estab- 
lished. It was developed in experimental programs 
conducted in the 1950s and early- 1960s by Grose 
and Taylor 1 , Schubauer and Spangenburg 2 , Spangler 
and Wells 3 , and Taylor 4 . Pearcy 5 summarizes much 
of the eady work conducted on vortex generator use. 
2. Secondary flow control. The rotational orientation 
of the vortex may also be used to counter boundary 
layer thickening due to the cumulative convective 
effects of secondary flows. This problem is typi- 
cally encountered in the ducted components of the 
aircraft propulsion system. Here the control mecha- 
nism originates with both the normal and transverse 
(left and right) components of the crossflow induced 
by the vortex. Recent experimental work in diffus- 
ing S -ducts 6, 7 have demonstrated the effectiveness 
of applying vortex generators to control secondary 
flows. 

Since the vortex generator is an obstruction pro- 
jecting into the fluid stream there is a drag penalty or 
loss of flow momentum resulting from its use. The goal 
of designing fluid components with vortex generators is 
thus one of maximizing component performance through 
one or both mechanisms of flow control described above, 
while minimizin g the number of devices used. This op- 
timization process presents practical difficulties due to 
the large number of parameters involved in a general 



Figure 1 Flow conditioning with vortex generators. 
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configuration of vortex generators. Consider Figure 1 
which illustrates an array of vane-type vortex generators 
in use on the wing surface of an aircraft First we must 
consider the geometric aspects of the individual vortex 
generators. These include chordlength, c, span, h, angle- 
of-attack with respect to the freestream direction, a, and 
other parameters related to the vane cross-section and 
axial profile. Secondly we also need to consider the ge- 
ometry of the array formation. This would include the 
spacing, s, between the individual elements in the array, 
the axial location of the array, zj, on the wing, and the 
rotational nature of the array, i.e. whether the vortices 
are shed counter or co- rotating with respect to imme- 
diate neighbors. Finally, the flow conditions impinging 
on the array (freestream Mach number, Reynolds num- 
ber, and boundary layer thickness for example) round out 
our list Limited experimental parametric studies have 
been conducted to develop effective strategies for vor- 
tex generator use in the subsonic diffuser flows described 
in References 6, 8, and 9. Of greater potential utility 
is the application of computational fluid dynamics to the 
optimization problem. In particular, fast component de- 
sign codes based on solutions to the parabolized (or ‘"re- 
duced'’) Navier-Stokes equations have the ability to run 
through large parametric ranges in a time and cost effec- 
tive manner. In the one notable example developed at 
NASA-Lewis, a code known as RNS3D (a modification 
of the PEPSIG RNS code based on the stream funcdon- 
vorticity formulation of the Navier-Stokes equations and 
developed for internal flows on nonorthogonal grids by 
Levy, Briley, and McDonald 10 ) has recently been used 
to optimize vortex generator use in a Defense Research 
Agency (DRA)-designed diffusing S-duct. Test data ob- 
tained at DRA verify the remarkable reductions in total 
pressure distortion predicted by the code. 7 

In these design codes vortex generators are not mod- 
elled as a part of the flow surface geometry of the duct 
but are represented by modelling the effects of the shed 
vortices on the flowfield. The mathematical models used 
are typically exponential equations representing the con- 
vective and wake properties induced by the vortex gen- 
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erator. These models are essentially boundary conditions 
on the velocity or vorticity fields applied at an axial loca- 
tion representative of the physical location of the vortex 
generators (typically at the trailing edge tip). This vortex 
generator modelling process consists of two distinct parts. 

In the first part a model of the convective and wake 
effects of the shed vortex or vortex array is constructed. 
To demonstrate how the model represents the convective 
nature of the vortex, consider Figure 2, which illustrates 
the secondary velocity field of an embedded vortex (data 
on the left, model on the right) shed from an airfoil- 
type vortex generator. This vortex generator is one of 
12 symmetrically placed vortex generators spanning the 
full inside circumference of a straight pipe. The data is 
acquired in the crossplane one chordlength downstream of 
the vortex generator trailing edges. The model flowfield 
is constructed from from a summation of 24 terms (12 
embedded vortices plus 12 image vortices to represent the 
convective influence of the wall) with each term having 
the form 11 : 


r t 


1 — exp 



( 1 ) 


where F{ is the measured circulation of the zth vortex or 
its corresponding image, wj nax is the measured peak or 
maximum value of the induced streamwise vorticity field, 
and Ri is the crossplane distance from the measured lo- 
cation of the vortex center. These quantities are known 
as “vortex descriptors” following the work of Westphal 
et al . l2 The model in Eq. (1) is based on the classical 
“ideal viscous” or Oseen vortex. 13 Other suitable models 
include potential vortices and Rankine vortices.* Exper- 
imental workers concerned with embedded vortices have 
often noticed the close similarity of observed crossplane 
vortex structure to these simple two-dimensional models. 
Potential vortices were first applied to results obtained 
in subsonic diffuser studies by Grose and Taylor. 1 Later, 
Pearcy 5 would develop embedded vortex interaction mod- 
els based on the inviscid representation of embedded vor- 
tex structure. Eibeck and Eaton 14 compared the structure 
of a single embedded vortex to that of a Rankine vortex. 
Studies by Pauley and Eaton 15 and Wendt et al . 11 * 13, 16 
have made comparisons to the Oseen model. To properly 
represent the flow losses induced by the use of vortex gen- 
erators a model similar to that described above is used 
to mimic the axial velocity deficit or wake downstream 
of the vortex generator. An example is the following 
relation developed by Kim and Patel 17 : 



( 2 ) 


Figure 2 A comparison of the transverse 
velocity field of an embedded vortex; data on 
the left, and a model constructed from the 
superposition of Oseen vortices on the right. 


The patchwork Rankine model uses the vortex vis- 
cous core radius, R C i , as a measure of vortex concentra- 
tion. Peak vorticity, circulation, and viscous core radius 
are simply related: R 7 ci — r</(irwj ,iax ). 
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Figure 3 The Internal Fluid Mechanics Facility of NASA Lewis. 


where v z is the primary or axial component of velocity, 
v cz is the ffeestream value of axial velocity. Av™ ax is the 
maximum axial velocity defect measured at the center of 
the viscous core of the vortex, and R ci is the previously 
noted core radius of the Rankine model. 

In the second part of the modelling process we relate 
the descriptors of the shed or modelled vortex to the 
collection of parameters illustrated in Figure 1, or in other 
words, define the following functional relationships: 

T = r(z,c,/i, a, profile, s, Af , ile, ...), 
u> max = u max (z,c,h y a, profile^, 6 , Af, Re , ...), 

R = R(z , c, h, a, profile, s,6,M y Re , ...), 

Av™ ax = A v™ ax (z y c, h y a, profile,*, 6, M y J!c, ...), 

(3) 

where “profile” refers to the shape or type of vortex 
generator, M is the ffeestream or core Mach number of 
the flow, and Re is a Reynolds number. 

In this experimental program the relationships in Eq. 
(3) are studied for one type of vortex generator (NACA 
0012). The vortex generators are mounted on the inside 
surface of a straight pipe, either in isolation or in counter- 
rotating arrays. Pipe core flow conditions are subsonic 
with a nominal turbulent boundary layer thickness to 
pipe diameter ratio, 6/d = 0.09. Three-dimensional 
mean flow velocity data are acquired in crossplane grids 
located at regular axial intervals downstream of the vortex 
generator or vortex generator array. This is done using 
a rake of seven-hole probes. 

Previous reports have considered the influence of 
vortex generator geometry and axial decay behavior. The 
initial vortex structure and decay of an array of counter- 
rotating vortices shed from a symmetric arrangement of 
vortex generators was described in Reference 11. The 
influence of vortex generator geometry and the axial 
decay of an isolated vortex was the focus of Reference 18. 
During the course of these studies it became apparent that 
the spacing between the vortex generators had a strong 
influence on the initial structure of the shed vortices. The 
primary objective of this paper is to consider the effect of 
vortex generator spacing on the shed vortex descriptors 


and subsequent axial decay. This is done by mounting 
the vortex generators at an s value which is half of that 
reported in Reference 11. The results reported here are 
compared against similar results in References 1 1 and 18. 

Facilities and Procedures 

This study was conducted in the Internal Fluid Me- 
chanics Facility (IFMF) of NASA Lewis. This subsonic 
facility, as it was configured for this test, is illustrated in 
Figure 3. Atmospheric air is supplied from the surround- 
ing test cell to a large settling plenum equipped with hon- 
eycomb flow straighteners and screens. Upon leaving the 
plenum the airstream is accelerated through an axisym- 
metric contraction and into the test section duct The test 
section duct consists of a straight circular pipe of inside 
diameter d = 20.4 cms. Driving this flow is a down- 
stream sub-atmospheric exhaust chamber which is contin- 
uously evacuated by central air supply compressors. The 
Mach number range in the test section is between 0.2 and 
0.8 with corresponding Reynolds numbers (based on test 
section diameter) between 0.95 and 3.80 million. Mass 
flow rates are between 3 and 7 kgs/sec. More information 
on the design and operation of the IFMF may be found 
in the report of Porro et ai [9 

Figure 4 is a cut-away sketch of the test section duct 
A short section of straight pipe (labelled “inlet pipe” in 
Figure 4) connects the exit of the facility plenum contrac- 
tion to the duct segment containing the vortex generator 
array. Static pressure taps located on the surface of the 
inlet pipe allow the nominal core Mach number in the 
test section to be set and monitored. 

The test section segment containing the vortex gen- 
erator array is referred to as the “vortex generator duct”. 
The inside surface of the vortex generator duct (and hence 
the attached vortex generator array) rotates about an axis 
coinciding with the test section centerline. 360 degrees 
of rotation is possible. The rotation is driven by a mo- 
tor and gear located in an airtight box above the vortex 
generator duct. 

The vortex generator is illustrated in Figure 5. These 
models consist of a NACA 0012 airfoil-shaped blade 
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Vortex Generator 



Figure 4 A cut-away sketch of the test section duct 


mounted perpendicular to the surface of a base plug. The 
surface of the base plug is contoured to the inside radius 
of the vortex generator duct. Both blade and plug are 
machined from an aluminum alloy. 

The cylindrical coordinate system used in this study 
originates in the vortex generator duct, z = 0 coincides 
with the trailing edges of the vortex generator array. 
The axial location of downstream (r, 0) crossplanes is 
given in terms of blade chord. For example, the first 
survey location possible is the (r, 0) crossplane located 
at z = 0.38c. 

The test section segment downstream of the vortex 
generator duct is non-rotating. This duct segment is la- 
belled “instrumentation duct” in Figure 4. The flowfield 
measurements are acquired in this duct through the use of 
a radially actuated rake-probe indicated in Figure 4. To 
acquire data in an (r, 9) crossplane, the rake-probe is first 
actuated over a radial segment extending from the duct 
wall to the duct centerline. The vortex generator duct 
and vortex generator array are then rotated an increment 
in circumferential position, A0, and the radial survey re- 


peated. In this fashion, pie-shaped pieces of the flowfield 
are examined. A narrow slot running the approximate 
length of the instrumentation duct allows the rake-probe 
to be located at various downstream crossplanes. A set 
of slot-sealing blocks determines the allowable axial lo- 
cation of survey crossplanes, 

z { - 0.38c, 1.00c, ..., 10.38c 
14.75c, 15.38c, ..., 24.75c. (4 ' 


The rake-probe consists of 4 seven-hole probe tips 
spaced 2.54 cms apart These probes are calibrated in 
accordance with the procedure outlined by Zilliac. 20 The 
flow angle range covered in calibration is ±60° in both 
pitch and yaw for the probe tip closest to the wall. The 
calibration range for the outer 3 tips is approximately 
±30°. Uncertainty in flow angle measurement is ±0.7° 
in either pitch or yaw, for flow angle magnitude below 
35 degrees (pitch and yaw flow angle magnitude did not 
exceed 35 degrees in this study). The corresponding 
uncertainty in velocity magnitude is approximately ±1% 
of the core velocity, v C2 . 



Top view. 

Figure 5 The airfoil vortex generator. 
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Experimental Results 



Figure 6 The vortex generator array. 

Figure 6 is a crossplane view of the duct (looking 
upstream) showing the array of vortex generators. The 
vortex generators are equally spaced circumferentially at 
mid-chord. A 9 b between the mid-chord position of ad- 
jacent blades is 15°. Each vortex in this array has a 
“domain-of-influence” in the crossplane. The circumfer- 
ential extent of this domain is defined by the behavior of 
the circumferential velocity component, v e , near the wall. 
For the interior vortex number 5 in Figure 6, for example, 
the flow near the wall and above the core of the vortex is 
moving clockwise {v$ negative). The flow above neigh- 
boring vortex core no. 6 is moving counter-clockwise 
(ve positive). Between these two vortex cores vg van- 
ishes along a roughly radial line. The location where 
this line intersects the wall is referred to as a “stagna- 
tion point” in vg. The circumferential extent of no. 5’s 
domain thus lies between the stagnation points in vg and 
is about equal to the spacing between vortex generators, 
i.e. about 15°, although this increases somewhat as the 
array develops downstream. The circumferential extent 
of the survey grid in this study is meant to capture the 
full domain-of-influence of one vortex thereby isolating 
the viscous core and allowing the vortex descriptors to be 
evaluated. The survey grid is between 20 and 35 degrees 
in extent and captures all of vortex no. 5 as shown by 
the shaded area in Figure 6. Note that a portion of vortex 
no. 6 (and, to a smaller extent, no. 4) is also captured. 

As discussed in References 11 and 18 crossplane 
grid resolution is dependent on the size of the viscous 
core. For axial locations upstream of the 2 = 6 chord 
crossplane. Ad = 1° and Ar = 1.3 mm is sufficient for 
an accurate assessment of circulation and peak vorticity. 
For 2 > 6 chords, Ad = 1.5° and Ar = 1.7 mm is 
sufficient 


Table 1 and Figures 7a-h summarize the results of 
this study. Table 1 lists the test conditions and derived 
vortex descriptors for each crossplane illustrated in Figure 
7. In the crossplane plots of Figure 7 the radial axis 
represents distance from the wall, in centimeters, and 
the circumferential axis represents angular position in 
degrees. The axial position of each crossplane is listed 
in units of vortex generator chordlength. The vortex 
descriptors (for vortex number 5) are listed below each 
set of crossplane plots in Figures 7a-h. 

The vortex descriptors originate from the transverse 
velocity data in the crossplane. The vector plots on the 
left-hand side of Figures 7a-h are the measured trans- 
verse velocity data. The transverse velocity data is first 
converted to streamwise vorticity data following the re- 


lation: 


6 vg vg 1 6v t 

Uz = 77 + 7 “ r 69 ' 


(5) 


where (u r ,w«) are the transverse components of velocity 
in the radial and circumferential coordinates, respectively. 
Finite difference formulas are used to represent the spatial 
derivatives in Eq. (5). The resulting streamwise vorticity 
fields are plotted on the right-hand side in Figure 7. Solid 
and dashed contour lines represent negative and positive 
vorticity, respectively. The contour increment for lines 
not labelled is ±3000 sec -1 . u mac is located at some 
grid point having coordinates (r e , 6 C ). The vortex circu- 
lation is calculated by first isolating the region of core 
vorticity in the data field. This is done by referring to 


A0 fc 

<r c ,e c ) 

z 

circulation 

peak vorticity 

deg 

cms, deg 

chord lengths 

m 2 /sec 

sec 1 

| Axial Deca 

y of the Narr 

1 

! 

>rtex Array] 


15 

(9.3,70) 

1.000 

-0.893 

-41752 

15 

(9.2*9) 

1*25 

-0.859 

-38894 

15 

(9.1*9) 

2350 

-0.830 

-36712 

15 

(8.7*9) 

3300 

-0.756 

-26869 

15 

(8.4*9) 

4.750 

-0.671 

-16763 

15 

(8.1,70) 

6.000 

-0.598 

-10161 

15 

(7.9,71) 

8300 

-0.457 

-3856 

15 

(7.6,71) 

10375 

-0338 

-2329 

| Initial Structure of the Widely-Spaced 

Vortex Array! 


30 

(9.3,92) 

1.000 

-0.722 

-35714 

| Initial Structure of the Isolated Vortex 



— 

(93,91) 

1.000 

-0.643 

-33189 

| Uncertainty] 





* (0.04,03) 

— 

t 0.009 

+ 300 


Notes: 

X. Descriptors for narrowly-spaced array are for vortex #5. 

2. See Figure 5 for the geometry of the vortex generator. 

3. Core velocity = 81 m/sec, boundary layer thickness (at z = 0) = 1.78 cms. 

4. Crossplane location of vortex generator trailing edge tip (for vortex 
generator shedding vortex #5): r Q - 9.2 cms, O 0 = 70 deg. 

Table 1 - A summary of test conditions and results. 
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the contour plots of voiticity in Figure 7. A path en- 
closing the region of core vorticity is defined. The outer 
boundary of the core is taken to be the location where 
strcamwise vorticity is 1% of w max . The circulation is 
then calculated according to: 

f v dl, (6) 

where V is the velocity vector in the crossplane and 1 
refers to the path coordinate. By using closed paths com- 
posed of line segments in the r or 9 coordinate direc- 
tions the circulation is easily determined. Un rtainty 


estimates for all vortex descriptors are given in Table 1. 
These are derived by combining the uncertainties in mea- 
sured velocities and probe location in accordance with the 
procedure outlined by Moffat 21 

Tj)c interaction between the vortex and the turbulent 
pipe boundary layer is depicted with contour plots of 
primary velocity ratio, v z /v cz . These are the middle 
plots in Figure 7. v cz is the core velocity of the pipe 
and is listed in Table 1. The crossplane profile of the 
vortex generator is indicated by the shading in the plot 
of streamwise vorticity on Figure 7h. The comer marked 
with the dashed cruciform is the trailing edge tip. 


a) z= 1.00c 



Center of Vortex: r = 93 cms, theta = 70 deg, Circulation = -0.89 m2/sec, Peak Vorticity = -41752 (1/sec) 


b) z= 1.62c 



Center of Vortex: r = 9,2 cms, theta = 69 deg, Circulation = -0.86 m2/sec. Peak Vorticity = -38894 (1/sec) 


c) z=2.25c 



Center of Vortex: r = 9.1 cms, theta = 69 deg, Circulation = -0.83 m2/sec, Peak Vortidty = -36712 (1/sec) 


Figure 7 Velocity and streamwise vorticity results. (Contim — ) 
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h) z=10 .38c 


M K *5 



Center of Vortex: r = 7.6 cms, theta = 71 deg, Circulation = -0.34 m2/sec, Peak Vorticity = -2329 (1/sec) 

Figure 7 Velocity and streamwise vorticity results. 


Vortex Structure at z — 1 Chord 

For reference purposes Table 1 contains vortex de- 
scriptors derived from two additional test cases not illus- 
trated here. The descriptors of a vortex shed from a single 
vortex generator (referred to below as the “isolated vor- 
tex” test case) are from results detailed in Reference 18. 
Also listed are descriptors of a vortex contained within 
an array of 12 counter-rotating vortices. These vortices 
were shed from vortex generators having an initial spac- 
ing s = A&t> = 30°. The vortex generators in this array 
thus spanned the interior circumference of the test section. 
This array formation is termed “widely-spaced”. The ini- 
tial structure and development of the widely-spaced array 
is detailed in Reference 11. 

The axial position at z = 1.000 chord is our refer- 
ence crossplane due to the observation (from the widely- 
spaced array results) that the axial profile of w max peaks 
here. Figure 7a illustrates the measured domain of vortex 
no. 5 at zjc = 1.000. Vortex 5 is in close proximity to 
its near-neighbor, vortex 6, and the two counter-rotating 
vortices form an “upflow pair”, i.e. the flow between 
vortices is directed away from the wall. The transverse 
secondary velocities are strongest in the radial coordinate 
in the region between the vortex cores. This differs from 
the structure observed for both isolated and widely-spaced 
array test cases where the strongest transverse velocity is 
circumferential and exists between the core and wall at 
zjc = 1.000. 

The circulation results listed in Table 1 show that 
the shed vortices strengthen as the spacing between vor- 
tex generators is decreased. We believe this strengthening 
results from localized areas of flow accelerating and de- 
celerating through and about the blockage created by the 
vortex generator array. A brief analysis is given below 
to support this statement 

It is also interesting to note that u max ~ T at 
z/c = 1.000 for all three spacing test cases. In fact, in 
all parametric variations considered in this program 11, 18 
(with the exception of vortex generator aspect ratio) cir- 
culation and peak vorticity are related in this manner. 


Figure 8 illustrates this for all descriptor data obtained 
at z/c = 1.000 and airfoil aspect ratio = 8h/irc — 2/tt. 
From a modelling standpoint, the as-of-yet unexplained 
relationship between w max and T for variation in aspect 
ratio is of critical importance. If it is possible to replace 
an initial z value of uj max with a relationship involv- 
ing T only, for example, the vortex generator modelling 
procedure (see Eq. 1) would be greatly simplified. 

Axial Decay 

Figures 9a-b illustrate the axial descriptor decay of 
the isolated vortex, the widely-spaced array, and the 
narrowly-spaced array vorticies. Figure 9a illustrates the 
decay of T, and Figure 9b the decay in u ? max . In general 
terms, circulation decay is gradual and has a linear profile 
while, in contrast, the axial decay of peak vorticity is 
rapid and exponential in profile. 

The amount of descriptor decay suffered by the 
narrowly-spaced array far outstrips the decay exhibited 
by the isolated vortex and widely-spaced array vorticies. 
An examination of Figures 7a-h suggests why this is so. 
Throughout its axial range of development the narrowly- 
spaced array vortex no. 5 maintains close proximity to 



Figure 8 Circulation versus peak vorticity. 
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its upflow partner. Contours of streamwise vorticity in- 
dicate that the vortex viscous cores are “in contact” for 
all crossplanes examined. For crossplanes downstream 
of z/c = 6.000 viscous core contact occurs between vor- 
tices in downflow pairs as well. Thus we suspect de- 
scriptor decay by vorticity diffusion to be the dominant 
decay mechanism. This decay mechanism has been pos- 
tulated in a previous study of tightly-spaced arrays of 
axial vorticies. 22 It was noted in Reference 22 that the 
decay rate of the descriptors was proportional to the radial 
gradient in streamwise vorticity at the crossplane location 
where the streamwise vorticity is zero (there will always 
be such a location between two counter-rotating vortex 
cores). The magnitude of this gradient increases as the 
distance between vortex cores decreases, i.e. as the ini- 
tial spacing s decreases. The present results confirm this 
model of vortex array decay. 


Modelling the Effects of Vortex Generator 
Spacing on Shed Vortex Circulation 

The effect of vortex generator proximity on shed 
vortex strength (or circulation) was first noted by Brown 
et al. % In Brown’s study, the array consisted of pairs of 



Figure 9 Axial decay of circulation and peak vorticity. 


vortex generators. Each pair was composed of two airfoil 
sections in close proximity, set to produce a downflow 
pair of embedded vortices. The spacing between neigh- 
boring pairs was much larger than the spacing between 
the vortex generators in each pair. The authors noted that 
when the distance between the vortex generators in each 
pair was decreased below a certain critical value the shed 
vortex strength decreased. This was termed “lift inter- 
ference”, but the effect was not quantified or analyzed. 
As discussed earlier, the results of this study indicate an 
opposite reaction; that is, counter-rotating vortices shed 
from an array of equally-spaced vortex generators are 
strengthened as the spacing decreases. We also note that 
the dominant vortex flow formation in this study consists 
of pairs of upflow vortices. 

As the spacing between vortex generators decreases, 
regions of local flow acceleration in the vicinity of the 
vortex generator array intensify. This assumption forms 
the basis for the following analysis. Let represent 
the circulation of a vortex shed from an isolated vortex 
generator. The circulation of a vortex shed from a vortex 
generator contained in an array of equal spacing Af?j is 

r a „: 

r a * =r iv +dT iv . (7) 

Dimensional analysis (using elements of inviscid airfoil 
theory) relates the circulation to the static pressure dif- 
ferential A P maintained about the airfoil: 


r iu 


cAP 
pVcz ’ 


( 8 ) 


where p is fluid density. Thus we write: 


dT iv = 


cdjAP) 

pVcz 


(9) 


Figure 10 diagrams a top view of an array of 4 equally 



Figure 10 A top view of the vortex generator array. 

spaced vortex generators similar to the formations consid- 
ered in this study. The airstream moves from bottom to 
top in Figure 10. The two middle vortex generators create 
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an upflow pair of vortices and the pair of vortex genera- 
tors on either side of this array produce downflow pairs. 
Counter to what we might expect intuitively, the static 
pressure rises between the two middle vortex generators. 
To see that this is so, consider what occurs when the 
spacing is decreased to the point where the trailing edges 
meet A region of stagnation pressure is established.* 
Correspondingly, static pressure decreases between the 
vortex generator pairs on either side of the array. Thus the 
effect of neighboring vortex generators is to promote an 
enhanced pressure differential d(AP) = A P a ~ P u — Pd 
about the chordlength of the vortex generators. This pres- 
sure differential is approximated inviscidly by assuming a 
mass flux balance between the leading and trailing edges 
of the array such that the axial velocity at the array trail- 
ing edge is (nominally) the core value, v cz . Thus: 


A P a = 



(/? + ) 4 -(/?-) 4 ] 
(0-) 2 (0 + f }' 


( 10 ) 


where: 


d 

/?_ = - A0j — csin a, 
p+ = -A0& + csina > 


( 11 ) 


and so: 


(fTfu — Kg 


VczC 

2 


O 3 _) 2 0?+) 2 j 


( 12 ) 


The resulting pressure differential is likely to be much 
less severe than that predicted by Eq. (10), due to 
the finite aspect ratio of die vortex generators and the 
resulting three-dimensional structure of the flowfield. The 
non-dimensional constant K a is introduced in Eq. (12) to 
account for this. Upon comparison to the circulation data 
collected in the three spacing test cases (at zjc — 1.000) 
we find K a & 0.0274. Following Eq. (12) we construct a 
plot of dTi v versus A0*. This is illustrated in Figure 11. 


Such a pressure field occurs in the vicinity of a vor- 
tex generator device known as a “Wheeler Wishbone”. 23 
Studies 24 * i3 indicate that a strong pair of upflow vortices 
are shed from such devices. 


dl/y m 2 /sec 



Figure 11 Circulation increase 
versus vortex generator spacing* 

Summary 

An experimental study is conducted to exam- 
ine the crossplane flow structure and axial decay of 
counter-rotating longitudinal vortices shed from an ar- 
ray of equally-spaced symmetric airfoil vortex genera- 
tors mounted on the surface of a straight pipe. The air- 
foils have a NACA 0012 cross-sectional profile and are 
mounted at an angle-of-attack a — ±16°. Impinging flow 
conditions are subsonic, incompressible, and the turbulent 
boundary layer thickness to pipe diameter ratio is approxi- 
mately 0.09. The vortex generator span to boundary layer 
thickness ratio is approximately 0.60. Vortex structure 
and decay are quantified using crossplane circulation and 
peak streamwise vorticity. The primary goal of this study 
is to examine the effect of spacing between the parent vor- 
tex generators. Crossplane velocity measurements and 
derived streamwise vorticity fields are presented for an 
array of vortices shed from vortex generators spaced 15 
degrees apart Initial structure (recorded at a crossplane 1 
chordlength downstream of the array trailing edges) and 
axial decay results are compared against similar results 
for an array of vortex generators spaced 30 degrees apart, 
and an isolated vortex shed from a single vortex genera- 
tor under identical impinging flow conditions. When the 
spacing between vortex generators is decreased the cir- 
culation and peak vorticity of the shed vortices increases. 
Analysis indicates this strengthening results from regions 
of fluid acceleration and deceleration in the vicinity of 
the vortex generator array. Decreased spacing between 
the counter-rotating vortices also produces increased ax- 
ial decay rates of circulation and peak vorticity. 
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